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ABSTRACT 

Recent full-sky maps of the Galaxy from the Fermi Gamma-Ray Space Telescope have revealed a 
diffuse component of emission towards the Galactic center and extending up to roughly ±50 degrees in 
latitude. This Fermi "haze" is the inverse Compton emission generated by the same electrons which 
generate the microwave synchrotron haze at WMAP wavelengths. The gamma-ray haze has two 
distinct characteristics: the spectrum is significantly harder than emission elsewhere in the Galaxy 
and the morphology is elongated in latitude with respect to longitude with an axis ratio ks2. If 
these electrons are generated through annihilations of dark matter particles in the Galactic halo, this 
morphology is difficult to realize with a standard spherical halo and isotropic cosmic-ray diffusion. 
However, we show that anisotropic diffusion along ordered magnetic field lines towards the center of 
the Galaxy coupled with a prolate dark matter halo can easily yield the required morphology without 
making unrealistic assumptions about diffusion parameters. Furthermore, a Sommerfeld enhancement 
to the self annihilation cross-section of ^30 yields a good fit to the morphology, amplitude, and 
spectrum of both the gamma-ray and microwave haze. The model is also consistent with local cosmic- 
ray measurements as well as CMB constraints. 
Subject headings: 



1. INTRODUCTION 

With the first year data release, the Fermi Gamma- 
Ray Space Telescope provided a wealth of new in- 
sights and detail of the gamma-ray sky. The en- 
ergy range and angular resolution of the Large Area 
Telescope (LAT) on board Fermi has significantly ad- 
vanced the understanding of many areas of gamma- 
ray astronomy, from point source studies like pulsars 
JAbdo et all l2009al l2010et ISaz ParkinsoneFaLf 1201(1 
and blazars JAbdo et al.l l2009bl , l2010dH3 IH). to diffuse 
emissions from the extragalactic gamma-ray background 
JAbdo et al.ll2010dlAckermann et aT1l2010U lAbdo et al " 



2010al) and the interstellar medium (ISM) JAbdo et al 



2009d iPorter et al.|[200l jStrong et atll2010h . 

Recently. iDobler et al.l (I2OIOD assembled full-sky maps 
of the Galaxy using the published raw photon data from 
Fermi from several hundred MeV up to several hundred 
GeV. These maps of gamma-ray emission from the dif- 
fuse ISM are produced primarily through three processes: 
cosmic-ray (CR) protons collide with the ISM producing 
7T° particles that decay to gammas, bremsstrahlung from 
CR electrons (and positrons) colliding with ions, and in- 
verse Compton (IC) scattering of starlight, infrared, and 
CMB photons by CR electrons. Because bremsstrahlung 
and 7T° emission are due to collisions of CRs with the 
ISM, these emissions are highly spatially correlated with 
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other maps of the int erstellar medium like the dust col- 
umn density map of ISchlegel et al.l J1998D . Since the 
IC emission is generated by interactions of CR electrons 
with the interstellar radiation field (ISRF) there is not a 
good morphological tracer of this emission at other ener- 
gies. However, CR electrons are primarily accelerated in 
supernova (SN) remnants and so their injection morphol- 
ogy should be very disk-like. Although diffusion effects 
are important, for isotropic diffusion through the Galaxy 
the resultant IC emission should also be very disk-like. 

Using template fitting techniques to morphologically 
regress out the emission from 7r°'s, bremsstrahlung, 
and IC from disk electrons from the Fermi maps, 
IDobler et all J2010I ) found an excess "haze" of IC emis- 
sion towards the Galactic center (GC) extending ±50 
degrees in latitude and with an axis ratio of roughly 
2.0. This Fermi haze is the gamma-ray counterpart 
to the microwave haze observed by the Wilkinson Mi- 
crowave A n isotrop y Pro be (WMAP) as desc r ibed i n 
iFinkbeinerl J2004al) and IDobler fc Finkbeinerl J2008aft . 
At WMAP wavelengths, the same electrons which gen- 
erate the Fermi IC haze interact with the Galactic mag- 
netic field to p roduce synchrotron microwaves. Recently 
Su et al. ( 2010) reconsidered the morphology, arguing for 
a "bubble" -like structure. Nonetheless, for reasons out- 
lined in $2] we use the "haze" moniker throughout this 
paper, although we are considering effectively the same 
gamma-ray signal. 

In both the gamma-ray and synchrotron cases, the 
haze emission is significantly harder than elsewhere in 
the Galaxy, implying that the electrons which produce 
the haze have a harder spectrum than the electrons ac- 
celerated and diffused through the Galactic disk. In fact, 
the required electron spectrum (number density per unit 
energy) is roughly dN/dE oc E^ 1 at high energies which 
is significantly harder than electrons generated by SN 
shock acceleration after taking into account diffusion ef- 
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fects. In that case, the steady state spectrum is closer to 
dN/dE cx E- 3 . 

The identification of the haze in both the WMAP and 
Fermi data imply that the haze is both real and that 
the underlying electron spectrum is very hard. It is 
this hard spectrum and the diffuse elongated morphol- 
ogy that are the defining characteristics of the emission, 
and any proposed origin for the electrons must match 
both of these features. For example, several authors 
have studied the connection between the haze electrons 
and young and middle aged p u lsars (IZhang et al.l [20091: 
Faucher- Giguer e fc Loebll20ia : iMcQuinn fe Zaldarriagal 
2010). The morphology however of the diffused electrons 
accelerated in pulsar winds would also be very disk-like 
and would not match the morphology 8 . Others have 
tried to reproduce the haze emission with a combina- 
tion of increased SN rate and modified diffusion parame- 
ters (IMcQuinn fc Zaldarriagall2010t iGebauer fc de Boerl 

2009) , but this also cannot produce the observed mor- 
phology or the observed spectrum, even including possi- 
ble reacceleration effects. Lastly, the re has been specula- 
tion that both the gamma-ray haze dLinden fc Profumol 

2010) and the microwave haze (jMertsch fc Sarkarll2010T ) 
are due to imperfect template subtraction, however nei- 
ther of these criticisms has been able to produce the 
morphology or the spectrum (amplitude and shape) of 
the observations using simulations. Furthermore, the 
gamma-ray haze is visible in the Fermi sky maps without 
performing any template fitting demonstrating that it is 
clearly a real structure. 

This work builds upon previous studies of the haze 
which explore the possibility that the haze electrons are 
generated throug h dark matter (DM ) annihilations in 
the Galactic halo. IFinkbeine r (2004b) originally showed 
that the microwave haze morphology and spectrum in 
the WMAP 1-year data was reasonably well matched by 
a DM model with a particle mass of M x ~ 100 GeV and 
with a self annihilation cross-section (av) ~ 3 x 10~ 26 
cm 3 /s which is roughly that required to yield the ob- 
served relic density of DM Q,dm ~ 0-23 if the DM parti- 
cle is a thermal relic of the Big Bang. 

However, initial data from Fermi of the inner Galaxy 
suggested that the IC emission from the haze elec- 
trons extended up to at least ~ 200 GeV implying a 
DM particle mass of closer to ~ 1 TeV. Since the an- 
nihilation rate is proportional to the number density 
squared, this requires a (av) roughly 100 times the ther- 
mal relic value in order to match the data. With light 
force carriers, a "boost fac tor" of 100 in the Galac- 
tic halo is easily obtai nable (|Arkani-Hamed et al.l 120091 : 
PosdcIov & Ritz 2009) via the Sommerfcld mechanism 
(ISommerfeldlll93lHHisano et alj|2005l 12004 ICirelli et al.l 
l2008HLattanzTfcS ilk 2009), in which (av) increases with 
decreasing relative velocity up to some satur ation value, 
while still producing the corre ct relic density (jFeng et al.1 
l2010t IFinkbeiner et al.l I2010D . Such a particle model 
is also consistent with local electron and positron 
CR anomalies observed by the Payload for Antimatter 
Exploration and L i ght-nuclei Astrophys i cs ( PAMELA ; 
iPicozza et all l2007t lAdriani et all 120091 |2010D satellite 

8 Millisecond pulsars in th e galactic halo may c ontribute to the 
haze signal at some level (see Malyshcv et al. 2010), but their mor- 
phology would also likely be spherical instead of significantly elon- 
gated in latitude. 



and Fermi (lAbdo et al.ll2009dt lAckermann et alj|2010al) 
as sho wn bv ICholis et al.l (|2009af ) and ICholis fc Weinerl 
( 2009). A model independent fit to all of the d ata (gam- 
mas, microwaves, and CRs) bv lLin et all (|2010l ) confirms 
that the injection spectrum must be E 2 dN/dE cx E 2 
which is broadly consistent with the spectrum of a Som- 
merfeld enhanced DM annihilation scenario in which the 
main products are leptons. 

These works have shown that the amplitude and spec- 
trum of the haze are easily reproduced with a DM parti- 
cle annihilation model; but here we are concerned pri- 
marily with the morphology. The morphology of the 
gamma-ray haze is the most difficult aspect to model 
since the haze is significantly elongated in latitude with 
respect to longitude. In fact, the geometry is impossible 
to realize with disk-like (or, as we show in 21 spherical) 
injection, ruling out SNe or pulsars as a possible source. 

Such a geometry is also inconsistent with a spher- 
ical DM halo and isotropic diffusion. However, it is 
very likely that neither of these assumptions is accu- 
rate. Generically, DM N-body simulations of Milky 
Way sized halos i mply prolate ha l os wi t h an axis ra- 
tio of roughly 2 (iDiemand et all 120081: iKuhlen et all 
l2008t ISpringel et al.ll2008D and observations of the spa- 
tial distribution of Milky Way satellites imply a prolate 
halo oriented perp endicular to the Galactic disk (e.g., 
IZentner et alJ l2005). In addition, the presence of any or- 
dered magnetic field lines towards the GC implies that 
the electrons will not diffuse isotropically as they follow 
the fields. In 52]we discuss the morphology of the haze in 
more detail, and in $3]we outline our anisotropic diffusion 
model which produces a DM IC halo that closely resem- 
bles the observed morphology. In fJH we compare our 
model to the data (both the morphology, amplitude, and 
spectrum of the haze emission) and in <J5] we summarize 
our conclusions. 

2. HAZE MORPHOLOGY 

Prior to the release of the g amma-ray data, th e mi- 
crowave haze was de s cribed by iFinkbeinerl (|2004al ) and 
iDobler fc Finkbeinerl (|2008aD as being centered on the 
GC, roughly spherical, and decreasing in amplitude ap- 
proximately as 1/r where r is the angular distance to 
the GC. However, such a microwave signal is limited by 
the extent of the B-ficld off the disk. The Fermi data 
on the other hand clearly show that the haze is in fact 
elongated in latitude b and extends to \b\ ~ 50 degrees. 
Despite the lower angular resolution and signal-to-noise, 
the gamma-ray data give a more complete picture of the 
location of the haze electrons. The reason for the dif- 
ferent morphologies is that the synchrotron amplitude 
is proportional to the magnetic field strength while the 
IC is proportional to the ISRF. Since the magnetic field 
falls off quickly with distance above the Galactic disk 
while the CMB amplitude is latitude independent, the 
microwave haze is confined to lower latitudes compared 
to the gamma-ray haze. 

The detailed morphology of the gamma-ray haze close 
to the Gala c tic p lane is difficult to determine. In 
IDobler et al.l (|2010l ). three methods of template fitting 
were used: 1) the actual Fermi data fro m 1.0-2.0 GeV 
was u sed as a full-sky template, 2) the Schleg el et al.1 
tfl998h (SFD) dust map was used alone , and 3) the 
SFD dust map, the Haslam 408 MHz map (jHaslam et al.l 
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Fig. 1. — Upper left: The haze residual using the 4 templates fit (disk, SFD, uniform, bubble) defined in |Su et al.| (2010T). The haze 
residual in this case is very pinched in the center and resembles two "bubbles" ; however, note the significant regions of over-subtraction near 
the disk which pinch the haze towards the center. Upper right: The haze residual using the 3 templates fit uniform, GALPROP). 

Although there is more noise, there is now very little disk over-subtraction and the haze looks much more like an "oval" . Bottom row: the 
same residuals with hand drawn contours over-plotted to highlight the morphological differences. 

sion elsewhere in the Galaxy. 

The advantage of method 3) is that the absolute spec- 
trum of the haze can be well determined since the haze 
structure is not in the external templates. However, be- 
cause of line of sight variations in the ISM and cosmic-ray 
proton density, there will be, for example, variations in 
the ratio of ir° gamma-ray emissivity to total dust col- 
umn density. Thus the dust column map will not be a 
perfect tracer of the gamma-ray map. This is especially 
true in the inner Galaxy (within about 30 deg of the 
GC) and has the potential to significantly effect the per- 



1982), and a bivariate Gaussian haze template were used. 
Method 2) was not particularly successful at fitting the 
full sky data and left significant disk-like residuals as 
well as the Fermi haze. Method 1) and 3) were much 
more successful but gave very different haze morpholo- 
gies at low latitudes (|6| < 30 deg). In particular, using 
method 1) gives a haze which is more oval shaped while 
method 3) gives a haze which is more h ourglass o r "bub- 
ble" shaped (see Figure HJ). Recently, ISu et all (|2010f ) 
explored the bubble morphology of method 3) in detail 
and argued that this morphology may be indicative of a 
significant event towards the GC (e.g., accretion onto the 
central black hole) in the past. However, before ascribing 
a physical mechanism to the generation of the haze elec- 
trons which is dependent upon the haze morphology, it 
is important to determine what that morphology is and 
why the two methods differ. 

Both methods 1) and 3) have associated problems. 
Since method 1) takes differences of Fermi data at dif- 
ferent energies, any haze that is present in the lower en- 
ergy data is subtracted off of the higher energy data so 
that the specific spectrum of the Fermi haze cannot be 
uniquely determined. In addition, since the Fermi maps 
have somewhat low signal to noise, subtracting one map 
from another (which adds the noise in weighted quadra- 
ture while removing the signal) yields difference maps 
that can be quite noisy. 

On the other hand, method 1) has the advantage that 
it does not rely on external templates (like the SFD dust 
map for example) and so automatically takes into ac- 
count systematics like line of density effects in the ISM. 
In other words, the lower energy Fermi maps are a bet- 
ter morphological tracer of the higher energy Fermi maps 
than external templates. The fact that the haze residual 
remains in the difference is a statement that this emis- 
sion has a significantly harder spectrum than the emis- 



ceived haze morphology. To illustrate this point, Figure[2] 
shows the Fermi data from 2.0-5.0 GeV with the lSu et all 
(f2010h model for IC emission and varying amounts of the 
SFD map subtracted. When the SFD coefficient is small, 
the 7T° gammas are clearly under-subtracted. However, 
as the coefficient is increased, a clear "X" shaped over- 
subtraction becomes visible. This structure defines the 
"bubble" shape of the haze in method 3), and may be the 
root of the discrepancy between the two morphologies. 
That is, if the haze were actually oval shaped, it may ap- 
pear more hourglass shaped after over-subtracting this 
"X". 

It is important to note that this "X" is not a feature 
in the SFD map (with the exception of the upper-right 
and possibly lower-right edges) but rather is being over- 
subtracted because the projected ir° to dust column ratio 
is lower in that shape. 

Furthermore, it is quite possible that the environmen- 
tal conditions towards the GC which give rise to this "X" 
in gammas, produce similar features in X-rays and mi- 
crowaves. For example, a heating source towards the cen- 
ter could heat the gas leading to enhanced, harder x-ray 
emission and such a variation in the environment would 
affect the estimate of column density to spinnin g dust 
emissivity used by iDobler fc Finkbeinerl (2008a. bj) and 
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Fig. 2. — The Fermi data at 2.0-5.0 GeV minus the disk IC model of ISu et al.l d20Tof) plus varying amplitudes times the SFD dust 
map as a tracer of 7r° emission. As the SFD amplitude is increased a clear "X" shape or over-subtraction emerges towards the Galactic 
center. This is due to lower w° to dust column ratio in that shape towards the bulge likely caused by line of sight density variations of the 
ISM and cosmic ray protons. This X-shaped over-subtraction can make the oval shaped haze (right hand panel of Figure [1} appear more 
"bubble" -shaped (left hand panel of Figure[TJ. 



iDobler et al.l (|2009f ) to remove the spinning dust com- 
ponent at microwaves. This would have the affect of 
making both the gamma-ray and microwave haze more 
hourglass shaped due to the same ISM physics which 
generates an edge in x-rays. Without speculating further 
what this "X" structure is, we note that there is signifi- 
cant evidence for X-shaped bulges in other galaxies, and 
recent evidence from the 2MASS survey that there exist 
red clu mp populations in the Milky Way that follow this 
feature (|Mc William fc Zoccali|[20ia) . 

In the context of comparing the gamma-ray haze spec- 
trum and morphology to a signal generated by inject- 
ing electrons via dark matter annihilations, the "bubble" 
morphology seems difficult to obtain (or at the very least, 
seems more indicative of a transient event in the GC). 
However, we show below that an oval shaped haze (and 
even an hourglass shaped haze) is possible with DM an- 
nihilation when considering anisotropic diffusion effects. 
Regardless, the underlying morphology of the gamma- 
ray haze at low latitudes is an unsettled issue. We choose 
to compare our results to the oval-shaped morphology 
and show that method 1) plus a dark matter contribu- 
tion to the IC emission with anisotropic diffusion effects 
is consistent with the data. 

3. DIFFUSION MODEL 

Since the basis for any anisotropic diffusion scenario is 
that electrons travel along ordered field lines, our diffu- 
sion model must first assume a geometry for the ordered 
component of the Galactic magnetic field. From there, 
this magnetic field can be related to specific diffusion pa- 
rameters which appear in the diffusion equation. All of 
our calculations are d one by modifying the CR propaga- 
tion code GALPROP (Strong fc Moskalenkoll 19981 120011: 
Moskalenko et al.ll200alPtuskin et al.ll2006tlStrong et al l 
2007) to include anisotropic effects. 



3.1. Galactic magnetic field model 

Our magnetic field model consists of two components: 
an irregular magnetic field B lrr and an ordered magnetic 
field -Bord- The former is parameterized as an exponential 
disk, 

= B e^- r)/ri ~^ /zi , (1) 

where r and z are the radial and vertical distances from 
the GC respectively, and -Bo is the local value of the 
irregular component (i.e., at r = Rq ~ 8.5 kpc, the GC- 
sun distance). The ordered field is assumed to have the 
form, 



-Bord = B\e 



-r/r 2 -\z\/ z 2 



(l + Ke- r / r3 -W/ z ^ , (2) 



where Bi(l + K) is the amplitude of the ordered field 
at the GC, which i s based on the 3D field model of 
lOrlando et all <f20Toh . 
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Fig. 3. — B-ficld amplitude profiles versus radial distance r (left) 
and versus distance above the galactic disk z (right). 



The parameters B , B%, K, 
that we use are set by hand 
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Note. — Magnetic field morphologies and parameters for the IC signals plotted in Figure [5] Our fiducial model is Model 1 which 
generates an IC signal that roughly matches the Fermi haze morphology (see Figure [BJ>. 



appropriate IC geometry and agree with measured 
values of the Galactic magnetic field at distances 
greater than ~ 1 kpc from the GC (I Jansson et al 



Sofue k Fuiimotol fl98l H an 
20011 : Han 2002t iTinvakov"^ 



Beck 



Qiaol 119941 
Tkachevl I2002t iSun et al 
2008t iBrown et alj 120071 : iBeckj l2009t Ijaffe et all 120101: 
Nishivama et al.l l2010l see Figure [3]). The parameters 
of our fiducial model (Model 1) are shown in Table [TJ 
These give a local value for the total magnetic field of 
5.4 /iG and an ordered-to-total amplitude ra tio of ~0-62 
which agrees well with measured values (see !Beck1l2009L 
and references therein). 9 

3.2. Anisotropic diffusion 

The propagation of CRs through the ISM is governed 
by the diffusion equation, 



~dt 



dQnp) 
dE 



(3) 



where tp is the number density per unit particle momen- 
tum of CRs at time t and position x, b is an energy loss 
coefficient (dominated by synchrotron and IC in the case 
of electron CRs), Q is a source term due to the injection 
of electrons by DM annihilations, and D is the diffusion 
constant. It is this last parameter which must be mod- 
ified for the case of anisotropic diffusion, and so we are 
concerned with the term above. 

We solve Equation [3] using GALPROP on a cylindrical 
grid so that, 



r or or oz Oz 



(4) 



Typically, isotropic diffusion is assumed so that D is 
not a function of x = (r, z). However in our case Eq. U] 
generalizes to: 



f(DfiP) = l^-( rDrr ^ + rD rz ^- 
r or or Oz 

d ,„ dth „ dtb. 
+ ^-(D zz -f + D zr -f), 
oz oz Or 



(5) 



9 These parameters do give a somewhat high value of 89 [iG 
for the total field at the very center, r = z = kpc. However, 
we note that not on ly is this in agreement with the estimates of 
ICrocker et"all lIMToT ) who pi ace a lower limit of 50 fiG in the inner 
400pc from necessary synchrotron cooling to avoid violating exist- 
ing diffuse 7-ray bounds, but also the very center is well outside 
our region of interest. Our mask of the Galactic plane extends up 
to \b\ = 5 deg or \z\ m 0.75 kpc. Inside this region, our choice of 
B-field has little impact on our results and our value at the center 
is only due to our specific parameterization of the field which likely 
does not extend in to arbitrarily small distances. 



where D rr , D zz , D rz and D zr are functions of x = (r, z). 
For details of the implementation of this anisotropy in 
the GALPROP code, see Appendix 

All that remains is to relate the diffusion tensor coef- 
ficients Dm Dzz an d D rz = D zr to the magnetic field 
model. iParkerl (|1965h describes the propagation of parti- 
cles along ordered field lines in the presence of an irregu- 
lar component, and in this case, the diffusion tensor can 
be written, 



Di 



Do 



n 2 



(6) 



where Do is the diffusion constant for the isotropic case, 
Sij is the delta function, fl; is the cyclotron frequency 
due to the field pointed along the i-direction (f2j oc Bi 
and = flf + n 2 ), and v is the characteristic frequency 
of deflections by the irregular component (y oc B m ). In 
our case, we assume for simplicity that the ordered field 
is oriented perpendicular to the Galactic plane, B r = 
and B z = £? rd, so that D rz = D zr = 0. In this case, the 
diffusion tensor becomes, 



Di 



D x 



(1 



o 



(7) 



where i? ra t is the ratio of the ordered to irregular field 
and we have used the fact that Cl/v oc B OT d/Bi„. Note 
that, in the limit of -B rd 0, D rr = D zz = Do, and in 
the limit of B m — > 0, D rr —> as desired. The form of 
this diffusion tensor implies that adding an ordered field 
suppresses diffusion perpendicular to that field. 

For the diffusion tensor coefficient, we assume Do oc 
E~ 0A3 . However, in contrast to most studies involving 
GALPROP, we incorporate the dependence o f -Pp on B tot 
as well. In particular following I Strong et al.1 (|2007f) . 



D oc 



Arr 

Btnt 



x r sy 



By, 



Bi 



(8) 



and because B depends on position, Do = Do(r,z). We 
set the normalization to be the locally measured value at 
roughly the locally measured magnetic field amplitude 
if the field were completely irregular, so that our final 
diffusion coefficient can be written as, 



D = 2.0 x cm 2 /s 



5 ^G 



E 



4.0 GeV 



-0.5 



where the normalization is fixed by fitting to the local 
CR measurements. 

Taken together, Equations [7] and |H] completely de- 
fine our anisotropic diffusion model and reduce to the 



6 



isotropic case when -B rd - ► and B; rr — > constant. For 
more details about the dependence of diffusion on the 
magnetic field, see Appendix [B] 

Lastly we note that, in all of our models, we use a box 
height Lbox = ±20 kpc. This is not directly compara- 
ble to the usual box heights (~ 4 kpc) discussed in the 
literature, because the "free escape" of electrons outside 
the Galactic disk is taken into account by the spatial de- 
pendence of the diffusion tensor. This is in agreement 
with findings by th e Fermi team regarding diffuse IC 
away from the GC (|PorterM201ClD and is in fact a more 
appropriate box size. This also alleviates the problem of 
"squashed" morphologies that are typical of smaller box 
heights when the CR density at the boundary is set to 



3.3. DM annihilation model 

In Equation [21 the source term Q is the rate of e + e" 
injection by DM annihilations and is given by 



n< \ l i \ dN 



P(r, z) 



(10) 



where dN/dE is the injection spectrum and p is the 
Gala ctic DM halo. We assume a prolate Einasto (|Einastol 
fl965h halo, 



p(r, z) oc exp 



a/2 



(11) 



with z c /r c = 2.0 , z c = 27 kpc, and a — 0.17 
(jMerritt et al.l l2005h . The overall normalization is set 
so that the local DM d ensity is p{R G ,0) = 0.4 GeV/cm 3 
(|Catena fc Ullid [2010h . 

The injection spectrum dN/dE is governed by the 
specific particle model. I n our case, we use XDM 
(jFinkbeiner fc Weinerll2007[ ) as our fiducial model, with 
M x — 1.2 TeV, an annihilation channel XX ~^ 4 >( t ) : 
<p —> e + e~, and with bran ching ratio 1 (hereafter, XDM 
e ± ; see lCholis et al.ll2009ctf bh. In this model, </> is a vector 
boson with < 2m^ that is the force carrier responsi- 
ble for the velocity dependent Sommerfeld enhance ment 
(jArkani-Hamed et all [2001 iPospelov fc Ritzll2009t ). We 
do not include specific dynamics for the host halo, but 
we do assume that the velocity dispersion (and hence the 
Sommerfeld enhancement or "boost factor") as well as 
substructure contribution is flat with radius. We define 
this boost factor BF as, 



BF = 



3 x 10" 26 cm 3 /s' 



(12) 



This model for the DM particle has E 2 dN/dE oc E 2 
as required by t he CR, microwave, and gamma-ray data 
dLin et al.ll2010h . 

3.4. Fitting Procedure 

We follow a similar pro c edure as that outlined in Ap- 
pendix C of lDobler et all ((2010). Specifically, we gener- 
ate a synthetic sky map, 

S(E) = A loE x + A sp x G(E) + U(E), (13) 

where A\ oE and j4 gp are the amplitudes of the Fermi 
0.5-1.0 GeV map and the GALPROP map at mean 



energy E = \/EqEx respectively and U{E) is a uni- 
form background, and convert to a synthetic counts map 
/i(E) — S(E) x (mask) x (exposure). We then minimize 
the log-likelihood, 



ln£ 



[ki ln/ij - m - ln(fci!)], 



(14) 



where ki is the map of observed counts at pixel i, over the 
parameters A\ oE and A gp . When comparing maps at dif- 
ferent energies, it is important to smooth the templates 
and data to a common beam full-width half-maximum 
(FWHM). All of our maps use 1.6 years of data, are 
smoothed to 2 degrees, and for t he £ n L ? map, we u se 
only "front" converting events (see lDobler et al]|2010l h 

4. RESULTS 

Figure @ shows the GALPROP IC map for E = 3.0 
GeV and for various assumptions about the dark halo 
prolateness and anisotropic diffusion. For the case of a 
spherical halo with isotropic diffusion (completely tan- 
gled magnetic field), the resultant IC signal is largely 
spherical. The same is true for our anisotropic model 
with a spherical halo, implying that diffusion effects alone 
cannot create the observed morphology. In fact pro- 
late halos lead to IC morphologies which very closely 
resemble the haze morphology. In detail, we find that 
the prolate halo with isotropic diffusion is overly concen- 
trated towards the center and that the best morpholog- 
ical match to the data comes from using a prolate halo 
with anisotropic diffusion. 

The detailed assumptions on the B-field morphology, 
and thus on the spatial dependence of the diffusion, can 
have a strong effect on the observed morphology of the 
IC emission as is shown in Figure [5] where we present the 
IC maps at 3 GeV for four distinctively different B or d as- 
sumptions from those of Equation [2j The specific mag- 
netic field model can lead to various IC morphologies 
from more uniform to more centrally concentrated and 
from more elliptical to more circular. In addition, for 
fields with a strong ordered component towards r = 
kpc, forked morphologies (due to increased synchrotron 
losses towards r = kpc) are found. Interestingly, for 
relatively modest changes to our magnetic field param- 
eters, we can also reproduce an hourglass shap e remi- 
niscent of the "bubble" shape in ISu et all (|2010l ) . Note 
that all models use an identical prolate dark matter halo; 
the variations in shape are due exclusively to magnetic 
field effects on the diffusion and relative energy losses to 
synchrotron and IC. 

This IC emission is the combination of electrons scat- 
tering CMB, IR, and starlight photons. Each of these 
ISRF components has a distinct morphology, and so the 
IC emission from each will also have a different mor- 
phology. In fact, since the starlight and IR photons are 
mostly confined to the plane, the high latitude IC emis- 
sion is due primarily to scattering of CMB photons. This 
is borne out in Figure [6] which shows the morphology of 
each of the IC components. The starlight and IR IC 
photons are concentrated much more towards the GC 
while the CMB IC photons extend to much higher lati- 
tudes. Furthermore, it is interesting to note the distinct 
"bubble" -like morphology of the CMB component. The 
implication here is that template fits lik e those used in 
iDobler et all ((20101 ) and lSu et all ((2010ft . which use ex- 
ternal templates that are concentrated towards the GC 
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Fig. 4. — GALPROP IC at 3 GeV due to production by DM annihilations with different assumptions about the halo shape and 
diffusion model: a spherical Einasto halo with isotropic diffusion (upper left), an axis ratio 2 prolate halo with isotropic diffusion (upper 
right), a spherical halo with anisotropic diffusion effects (lower left), and a prolate halo with anisotropic diffusion effects (lower right). 
All plots are arbitrarily normalized to the same intensity at (£, b) = (0, 50) degrees. The spherical halos are clearly inconsistent with the 
haze morphology (see the right hand panels of Figure [TJ while the prolate halos provide a significantly improved fit. In particular, the 
anisotropic diffusion case gives a morphology that has both the observed axis ratio and concentration. 
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Fig. 5. — The same as the bottom right panel of Figure [4] but for several different models of the ordered magnetic field (see Table [T}- 
Different magnetic field models can lead to various IC morphologies including forked (top left, due to increased synchrotron losses towards 
r = kpc), circular and centrally concentrated (top right), circular and more uniform (lower left), and also more hourglass-shaped (lower 
right). See j4]for a description. 



could potentially absorb the starlight and IR morpholo- 
gies, while leaving the CMB morphology which appears 
more bubble-like. That is, if the intrinsic haze morphol- 
ogy is more oval-shaped, pulling out only the CMB com- 
ponent would may leave a bubble morphology. 



In Figure [7] we show the residual "haze" map, 



Ul\ = Ell - Sf 1 



+ A gp xG(E), 



E E E 

as well as the residual map, 

Eo Eq Eo 

As shown in the hgure, the three component model pro- 
vides a remarkably good ht to the data. There is some 
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Fig. 6. — The full anisotropic, prolate model at 3 GeV (top left) as well its individual components broken down into photons from IC 
scattered CMB (top right), IR (bottom left), and starlight (lower right) photons. All maps use the same stretch, normalization, and contour 
intervals. The CMB component in particular has a distinctly "bubble"-like morphology. Thus, when using template regression techniques 
to assess the underlying morphology of the haze, care must be taken not to regress out emission from IR and starlight components while 
leaving only the CMB component. 



residual over-subtraction due to the fact that the Fermi 
haze appears to have an "edge" at roughly \b\ ~ 50 deg. 
This feature cannot be reproduced exactly by our mod- 
els which tend to be slightly more diffuse. This lack of 
an edge pushes the fit to slightly over subtract the GAL- 
PROP haze c ontribution. Indeed, astrophysical models 
such as winds (|Crocker et alJl201ll ) or jets would also ei- 
ther not have an edge or, in the case of jets, likely have a 
shock heated edge with a harder spectru m which is not 
clearly seen in the data (|Su et al.l I20 10) . Despite this, 
our fit removes 96%, 89%, and 69% of the variance over 
pixels with |6| > 5°at E = 2-5, 5-10, and 10-20 GeV 
respectively (see Figure [7|). 

Lastly, we compare the spectrum of the observed Fermi 
haze to that produced by the IC emission from e gen- 
erated by the XDM electrons annihilation channel. We 
plot the H(E) emission in the window defined by \£\ < 20 
deg and 10 < |6| < 50 deg. This region is dominated by 
the Fermi haze and is relatively free of other foregrounds. 
When comparing the spectra, it is important to keep in 
mind that H(E) has the Fermi Eq[® map times A\ e(E) 
removed, and so in Figure [8J we show the intensity versus 
E for H(E) and G{E) - A loE {E) x Gj;g. 

Performing an independent fit of the Fermi and 
WMAP haze profile (intensity as a function of latitude 
south of the GC) we find that the required BF for the 
Fermi haze at 4 GeV is BF=24 while at WMAP 23 GHz 
it is a nearly identical BF=27, as shown in the bottom 
panels of Figure [51 In our calculations of the synchrotron 
radiation emission, we take into account the presence of 
both the ordered and the irregular B-field components. 
In the upper right panel of Figure [5J the predicted DM 
spectrum is plotted over the Fermi data assuming a BF 
of ^24. It is clear from the figure that the DM spectrum 
with this BF provides an excellent agreement with the 
data, especially taking into account uncertainties in the 



optical and IR ISRF at latitudes far above the plane. 
While the cross section in the inner galaxy is roughly 
a factor of three lower than that needed to explain lo- 
cal cosmic ray excesses, this could nat urally arise from 
a rad ius dependent velocity dispersion (jCholis fc Weinerl 
120091 ) . or from a depletion of substructure in the inner 
galaxy (jSlatver et al.ll201lD . 

5. CONCLUSIONS 

We have developed a model of Galactic cosmic-ray 
diffusion that incorporates both an ordered and turbu- 
lent magnetic field component. The ordered component 
results in anisotropic diffusion of cosmic-ray electrons 
along field lines. Combining this model of diffusion with 
dark matter annihilations in a prolate Galactic dark halo 
produces an inverse Compton gamma-ray signal that 
matches the morphology and spectrum of the observed 
Fermi gamma-ray haze. Namely, an oval-shaped haze 
with axis ratio « 2.0, extending up to |6| ~ 50 deg, and 
with a cosmic-ray injection spectrum E 2 dN/dE oc E 2 . 

The detailed morphology of the haze at low latitudes 
is still uncertain. We have shown that the dust-column 
to 7T° gamma-ray ratio is higher in an "X" shaped mor- 
phology towards the center of the Galaxy and that using 
a map of dust column like the SFD dust map as a tracer 
of 7T° gammas results in an over-subtraction of the "X" . 
The end result is that an oval-shaped haze may then ap- 
pear more "hourglass" or "bubble" shaped. Using the 
0.5-1.0 GeV Fermi map itself (which contains very little 
of the gamma-ray haze) as a tracer of disk emission at 
higher energies is immune to these line of sight effects and 
produces a more oval-shaped haze at the cost of noisier 
residuals. 

Regardless, a three component model of anisotropic 
diffusion with dark matter annihilations in a prolate halo 
plus the Fermi 0.5-1.0 GeV map plus a uniform back- 
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Fig. 7. — The haze (left column) and residuals (right column) of our three template fit in energy bins from 1.0 to 20.0 GeV. The haze 
maps clearly show the strong haze residual with an axis ratio Ri2 that is in reasonably good morphological agreement with our anisotropic 
dark matter model (see Figure \E\ top left panel). This is borne out in the residual maps which show a residual consistent with noise at 
latitudes above 20°. Close inspection reveals a slight over-subtraction towards the center due to the fact that our model does not explicitly 
include an "edge" at 6 f» ±50° as is seen in the data. 



ground provides an excellent fit to the data from 1-20 
GeV. The self-annihilation cross section required for the 
dark matter generated IC component is ~ 9 x 1CP 25 
cm 3 /s (boost factor ~ 30), which is easily obtainable 
via the Sommerfeld enhancement in our models and also 
produces the microwave haze. Furthermore, this boost 
factor is well within the bounds of thermal relic and CMB 
constraints (|Slatver et all 120091: IZavala et al.ll2010D . 

The most significant outstanding issues are the sharp 
"edges" of the haze at high latitudes and also the mor- 
phology of the haze at low latitudes. Sharp edges are not 
particularly expected with either a dark matter annihila- 



tion or astrophysical (such as winds or jets) mechanism, 
unless the spectrum at the edge is significantly hardened 
as does not appear to be the case. Magnetic confinement 
could potentially help both explanations, though care 
must be taken not to significantly synchrotron brighten 
the edges which are not seen in the WMAP microwave 
data. The low latitude morphology of the haze ("oval" 
versus "bubble" shape) may become more clear as more 
data are collected by Fermi. In particular, at high en- 
ergies, the disk fades much more quickly than the haze 
because of the softer spectrum of the disk, and so the 
low latitude haze may be revealed at high energies with 
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Fig. 8. — Upper left: A zoom in of the haze from Figure [TJ with the two black boxes indicating the area used to plot the spectrum 
shown in the upper right. The light and dark bands represent the spectrum of emission in the north and south boxes respectively. The 
spectrum from DM annihilations (solid lines) is quite consistent with the data especially when taking into account uncertainties in the IR 
and starlight components of the ISRF at high latitudes. The dashed and dot-dashed lines are for the same model with varying IR and 
starlight intensities (20% and 50% reduction in IR and starlight intensities respectively). The required boost factor is nearly identical to 
that found by fitting the radial profile south of the GC to the Fermi 4 GeV data (bottom left) as well independently fitting the microwave 
haze profile at WMAP K-band (bottom right). 



5 to 10 more years of data. 

Of course, there is the possibility of a hybrid scenario in 
which some event evacuates a cavity towards the Galactic 
center that is filled with high energy electrons from dark 
matter a nnihilation that are trapped by magnetic con- 
finement. ISu et al.l (|2010f) discount this possibility under 
the assumption that the dark matter signal would be 
more spherical, but we have shown here that this is not 
the case in general for triaxial halos. Injection from dark 
matter annihilation would also have the advantage that 
the hard spectrum can be obtained (as we have shown in 



this paper) and the injection is extended. Nevertheless, 
inside the edge, the haze appears to have a profile that is 
roughly flat in latitude above |6| > 30°. Such a projected 
profile seems nearly impossible to realize (either with as- 
trophysical or dark matter models) unless electrons pile 
up on the edges, though the naive expectation would be 
that the gammas would be limb-brightened which is not 
observed. If these features persist in future data, such 
hybrid scenarios are inevitable. 

Lastly, we point out that introducing a significant or- 
dered field as we have could potentially produce a sig- 
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nificantly polarized microwave signal. By design, our 
model docs reproduce the observed microwave haze in 
total temperature; comparison with the WMAP polar- 
ization data will be the subject of future work. 
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APPENDIX 

A. ANISOTROPIC DIFFUSION IN THE GALPROP CODE 

In GALPROP, the diffusion equation is solved through the Crank- Nicholson implicit method (Str ong et al.l [20071 : 
iPress et al.|[l992f ): 

~dt~ At ~ At +Ql > (A1) 

where i is the index of position (r or z) or momentum and ai^a/At are the Crank-Nicholson coefficients. In the case 
where D is homogeneous in space, these coefficients are: 

El = n 2ri Z Ar El = n 2n El = n 2ri + Ar (AD 

At 2r,(Ar) 2 ' At r,(Ar) 2 ' At 2r l (Ar) 2 ' 

for diffusion along r and 

a% D a 2 2D a 3 D 



At (Az) 2 ' At (Az) 2 ' At (Az) 2 



(A3) 



for diffusion along z. With the new terms from anisotropic diffusion (see Eq. [5]) the Crank- Nicholson coefficients 
become: 

Oi _ 2r. t - Ar D rri+1 - Drr^ a2_ _ „ 2r { 03 _ 2r a + Ar D r r l+1 ~ Axj-i 
At~ m 2n(Ar) 2 4(Ar) 2 ' At ~ m n(Ar) 2 ' At ~ m 2r<(Ar) 2 + 4(Ar) 2 ' [ ' 

for diffusion along r and 

Oi _ Dzz, Dzzj+i ~ Dzz z -! <X2 _ 2D ZZz a 3 _ D ZZi D Z z 1+1 - £W_i 
At ~ (Az) 2 4(Az) 2 ' At ~ (Az) 2 ' At ~ (Az) 2 4(Az) 2 

for diffusion along z, and where, as in the main text, we have taken D rz = D zr = (since B r = 0) for simplicity. 
When iterating recursively for a steady-state ip, GALPROP utilizes the fact that the r and z directions are separable, 
whereas if we take D rz = D zr 7^ 0, these directions are not separable. So while we use Eqs. IA4I and IA5I in practice, 
the general quantization of Eq. [5] is (superscripts are spatial indices for clarity) 



(A5) 



rr 1 rr i r r : '- r\ X 



r l 2Ar 2Az / 2Ar 

^rz 1 rz rz 1 zz zz \ r r 



(A6) 



-* 2Ar 2Az / 2Az 

. ^+i,j + ^«-i,j_ 2 ^.j ^J+i + y^i-i - 2^-? 

rV X (Ar) 2 ^ X (Az) 2 + 

Dli x 2a7az • 

B. DIFFUSION DEPENDENCE ON MAGNETIC FIELD 

Let us consider the generic case of an electron traveling in a magnetic field with both an irregular B m and ordered -B rd 
component. As the electron spirals around the ordered field lines with cyclotron frequency f2, there is a characteristic 
frequency v at which the electron is scattered from its path by the irregular component. In the case of a strong 
ordered component B ov< \ £>i rr , il >• v while for B m & <C B m , <C v. In other words, for strong ordered fields, the 
electron spirals around the field line many times before it is de flected by the irregular component. As noted in ^3.21 
this behavior is written in the diffusion tensor as (|Parkerlll965j ): 

/ v 2 Sa + fLfL \ , , 

ij = D ° ( J 2 + o 2 ) • ( B1 ) 

which, for the case of an ordered field completely along the z-direction, leads to the relation 

ElL oc (B2) 

Dzz l + B 2 ld /B? T 1 ; 

We wish to motivate Equation IB2I from the perspective of diffusion lengths A since that it the most direct measure 
of the diffusion coefficient, Dij — Ay-c/3 (where c is the speed of light). In our scenario then, 

hi „ r ^EH „ N: (B3) 
A rr ^gyr 

where r gyr is the gyroradius and N is the number of scatterings of the particle by angle <j> ~ B j TT /B to t, which is the 
inclination angle of the field lines from the direction of the mean field due to irregularities (see|Longair 20021). Note 
that, in the case of B m 3> i? rd, 4> is large and particles are deflected significantly from their initial direction within 
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one gyroradius while for B m <C i? rd, ^ C 1 and particles follow the field lines of the local ordered field. For the 
particle to scatter by ~1 radian, we need yN(j) ~ 1 which implies 



_££ r~j (f> 2 ~ B 2 ot /Bf m ., 



and since B tot = yjB* rd + Bf rl , 



as desired. 



^ ~ r 2 ~ a + bijbd 



C. LOW FREQUENCY RADIO EMISSION 
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Fig. 9. — Upper left: The full sky Haslam 408 MHz map. Upper right: the latitudinal profile of Haslam and our 408 MHz synchrotron 
GALPROP map for primary electrons injected via SN shocks using our anisotropic diffusion model (both with and without the ordered 
field component). Middle left: Haslam minus our 408 MHz GALPROP map. While the agreement is in general quite good, the map is 
slightly too concentrated towards the center. This is alleviated by removing the ordered field {middle right) indicating that more complex 
ordered field morphologies than our simple parameterization allows can provide a good fit to the data. In addition to being a good fit to 
Haslam, the model without an ordered field is a reasonable match to the Fermi haze, though it is a bit centrally concentrated {bottom 
panels). 



Do minated by synchrotron emission from electrons with energies ~few GeV, the Haslam 408 MHz map (Haslam ct al. 
1982) provides an excellent constraint on both the injection morphology of the primary electrons via SN shock accel- 
eration and also the magnetic field morphology. Thus, it is important to check that our anisotropic diffusion model 
maintains the "disk-like" shape of synchrotron from primary electrons at 408 MHz. Figure |9] shows the Haslam map, 
and the Hasl am map and the predicted 408 MHz synchrotron emission due to SN injection using our full anisotropic 
model. As in lLin et al.l (|2010l ). we normalize the model to Haslam by setting the total emission in the region \t\ < 10° 
and -90° < b < -5° equal. 

As shown in the radial profile panel (upper right), the agreement is very good, though the map difference (middle left 
panel) indicates that the stronger magnetic field in the center (due to the ordered component) may make the emission 
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somewhat steeper in the region ~l-2 kpc. However, we point out that, not only is the gamma-ray signal dominated 
by emission at higher latitudes, but we have used a very simple parameterization for the ordered field and additional 
field parameters (to lower the field within the inner 2 kpc) can remove the discrepancy. In fact, we can remove the 
ordered component altogether and the model comes into close agreement with both Haslam (middle right panel) and 
the Fermi data (bottom panels; also cf. Figure |4j. 



